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We study the charge correlations between charm mesons produced in 500 GeV 
7r~ —N interactions and the charged pions produced closest to them in phase space. 
With 110,000 fully reconstructed D mesons from experiment E791 at Fermilab, the 
correlations are studied as functions of the Dn — D mass difference and of Feynman 
x. We observe significant correlations which appear to originate from a combination 
of sources including fragmentation dynamics, resonant decays, and charge of the 
beam. 



In the production of a D meson the charge of pions formed nearby in 
phase space may be highly correlated with the charge of the D meson. For 
example, this associated pion can be produced in the fragmentation dynam- 
icsO of a c quark combining with a d from a d — d vacuum pair, forming a D + , 
and the remaining d forming an associated pion nearby in phase space. Thus, 
D + (D* + )tt~ , D~ (D*~)tt + , and D°tt + D°ir~ would be favored over combina- 
tions with opposite pion charges. Resonant decaystl and the charge of the beam 
particle are also sources of charge correlations. Using a n~ beam can lead to 
the association of both charm mesons and anticharm mesons with negative 
pions, especially in the forward (beam) direction. By comparing the charge 
correlations of different species of charm mesons and antimesons with associ- 
ated pions, and by studying them as functions of Feynman x (xp), one can 
hope to disentangle some of these processes. 

We use well identified D° ->• K~Tr + , D+ ->• K-n+ir+, and D*+ -> D tt+ 
signals (and their charge conjugate decays) from experiment E791 at Fermilab 
for this study. A silicon vertex detector and two Cerenkov counters are used to 
locate decay vertices and identify pions and kaons for this..analysis. Details of 
the spectrometer and event selection are found elsewhere.tj The final D meson 
sample sizes used for this analysis are 22587 ± 210, 24237 ± 216, 24569 ± 204, 
29649 ±238, 4997 ±84, and 6048 ±93 events, for D°, D°, D+ , D~, D*+, and 
D*~ , respectively. 

For each D found in an event, all identified pion tracks originating from 
the primary vertex are combined with the D. Among these combinations, the 
pion that forms the smallest invariant mass difference (Am m i n ) with the D is 
selected as the candidate associated pion. 
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Table 1: The xp- and Am-integrated correlation parameters a defined in Eq. (1) for the 
background-subtracted signals prior to correction, for the corresponding backgrounds, and 
for the signals after correction using the matrix correction technique. 



Charm 


Signal a 


Background a 


Corrected Signal a 


D° 


0.13 ±0.01 


-0.04 ±0.01 


0.12 ±0.03 ±0.04 


W 


0.18 ±0.01 


0.04 ±0.01 


0.42 ±0.02 ±0.03 


D+ 


0.18 ±0.01 


0.10 ±0.01 


0.45 ±0.03 ±0.03 


D~ 


0.08 ±0.01 


0.02 ±0.01 


0.03 ±0.03 ±0.04 


D*+ 


0.15 ±0.02 


0.08 ±0.03 


0.33 ±0.05 ±0.03 


D*~ 


0.08 ±0.02 


0.02 ±0.03 


0.15 ±0.05 ±0.04 



We define the correlation parameter a as 

n(m= T,Ni(D* v )-ZNi(D*- q ) m 
1 ; ~ ENiiDn^ + JZNiiDir-i) [ 1 

where q = +, -, -, +, -, + for D = D°, W , D+, D~, D*+, and D*~ , 
respectively, and ^ Ni(Dn q ) denotes the number of charm mesons for which 
the selected pion has the charge q. 

The correlation parameters for background-subtracted signals and back- 
ground regions are listed in Table 1. The signal correlations differ significantly 
from the background correlations. These raw correlations must be corrected 
for the occurrence of falsely reconstructed tracks (ghosts), momentum smear- 
ing in bins of Am m i n , detector acceptance, and dilution of the correlation by 
wrong track selection. We use a Monte Carlo simulation-of the experiment and 
the LUND event generator (PYTHIA 5.7/JETSET 7.3)14 to model these effects 
and extract the final results (Table 1), within a matrix formalism explained 
elsewhereta A fundamentally different model of hadron production could change 
the corrections; an effect not accounted for in Table 1. 

To investigate beam-related effects in more detail, we study the xf de- 
pendence of the D + and D* + correlations. In Fig. 1, we plot a as a function 
oi xf, for both particle and antiparticle for D + and D* + . We observe that 
a(D + ) rises slightly with xf but a(D~) falls sharply to negative values for 
xf > 0.2. In both cases, the D's are more likely to be associated with 7r~'s at 
high xp where beam effects seem to be important. There appears to be less 
dependence of a on if for the D*^. 

Detailed Monte Carlo studies suggest minimal or no beam-related effects 
when the combined particle and antiparticle correlations are computed. We 
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calculate the combined correlation parameters to be a(D°, D°) = 0.29 ± 0.02 ± 
0.03, a{D+,D~) = 0.21 ± 0.02 ± 0.03, and a(D*+,D*-) = 0.23 ± 0.04 ± 
0.03. These results indicate that fragmentation dynamics and resonant decays 
produce substantial correlations between D mesons and their associated pions. 
All three combined correlation levels are approximately equal. 

In summary, we observe significant production correlations between D 
mesons and their associated pions. Beam-related and other fragmentation 
effects play an important role. Further studies are needed to unravel the con- 
tributions. 
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Figure 1: The corrected correlation parameter a as a function of xp for (a)D+ and (b)D*+. 
The parameter a is defined in Eq. (1) in the text. The error bars correspond to the statistical 
and systematic uncertainties added in quadrature. 
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